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Abstract 

Abstract: High temperature expansions for the susceptibihty and the sec- 
ond correlation moment of the classical A^-vector model {0{N) symmetric 
Heisenberg model) on the sc and the bcc lattices are extended to order 
for arbitrary A^. For A^ = 2,3,4.. we present revised estimates of the critical 

parameters from the newly computed coefficients. 
PACS numbers: 05.50+q, 11.15.Ha, 64.60.Cn, 75.10.Hk 
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There has been a resurgence of interest in series expansions for the Statistical Mechanics 
of lattice models witnessed by the recent publication of several new remarkably long high- 
temperature (HT) and low-temperature (LT) series, in particular for the A^— vector model 

with = (the self-avoiding walk model) 0, with = 1 (the Ising model) 0] and 
with N = 2 (the XY model) [Q in 2- and 3-space dimensions. The best results, however, are 
still restricted to the = and 1 cases, where series are obtained by counting techniques 
which achieve maximal efficiency in low dimensions and only with discrete site variable 
models. Presently, on the sc lattice, the zero field susceptibility xi^i P) and the second 
correlation moment fi2{N] (3) are known to 0{(3'^'^) and 0(/3^^) 0,^ respectively for A^ = 0, 
to 0(/5i9) and 0(/5i5) respectively for A^ = 1, to 0(/3i^) |] for A^ = 2 and to 0(/?i^) 
PJ^ for any other A^. On the bcc lattice, x ^md have been computed to 0{l3^^) for 
AT = [|, to 0(/52i) for A^ = 1 [|ig,|n|], to 0(/3i2) for A^ = 2 m and to 0{13^^) for A^ = 3 



|r3| . Apart from the interest of an increasingly precise direct determination of the critical 
properties of the lattice models, there is no lack of other good reasons to undertake such 
a laborious calculation as a long series expansion: they include more accurate tests of the 
validity both of the assumption of universality, on which the Renormalization Group (RG) 
approach to critical phenomena is based, and of the approximation procedures required to 
produce estimates of universal quantities by field theory methods. In fact, for want of more 



rigorous arguments, as stressed in Ref . , a crucial test of the validity of Borel resummed 



e— expansions or fixed dimension (7— expansions |[T^-|17| is still provided by the comparison 
with experimental or numerical data. 

Here we present a brief analysis of a new extension from 0(/?^^) to 0(/9^^) of the high 
temperature (HT) expansions in zero field for the susceptibility and the second correlation 
moment of the A^-vector model both on the sc lattice and on the bcc lattice. More results 
to 0(/5^^) in d = 2, 3, 4, .. space dimensions, and a study of the second field derivative of the 
susceptibility x^^\N-,(3) = d^x/dH^ will appear elsewhere. We have determined the series 
coefficients as explicit functions of the spin dimensionality A^ by using the vertex renor- 
malized linked cluster expansion (LCE) method UTSf. Our calculation extends substantially 



previous work in that it provides coefficient tables of considerable length irrespective of the 
spin dimensionality and of the lattice structure: HT series for the general A^-vector model 
were previously available only for the (hyper)sc lattice in ci = 2, 3, 4 dimensions up to 0(/5^^) 

Concerning the LCE technique, we have found the following works particularly useful: 
the review papers [|1|], the = 1 computations |ll0|,|n],|l9|-|ll] and the more recent work by 



M. Luescher and P.Weisz (LW) [|], devoted to the model, the O(A^) symmetric P{(p^) lattice 
field theory, described by the partition function 

Z = Udn{ipf)exp[/3 0i ■ ^j]^ (1) 

where 0i is a A^-component vector. With the choice diji{0l) = 6{0j — l)dipi of the 
single spin measure, (||) reduces to the partition function of the A^-vector model, but also 
a broad class of other models of interest in Statistical Mechanics can all be represented in 
this form. LW have devised or simplified some algorithms required for the calculations, 
and have tabulated HT expansions of A^2; X^^^ the (hyper)sc lattices for the A^-vector 
model to 0(/?^^) ||§J^. Also starting from (|T]), we have extended the calculation to the class 
of bipartite lattices, in particular to the (hyper)sc and (hyper)bcc lattices. By redesigning 
the algorithms in order to take full advantage of the structural properties of the bipartite 
lattices and by writing an entirely new optimized code, we have significantly reduced the 
growth of the complexity with the order of expansion. Thus we have been able to push 
our calculation well beyond 0(/3^^) where LW had to give up. We can give a rough idea 
of the size and the complexity of the calculation by mentioning that over 2 • 10^ graphs 
enter into the evaluation of x and /i2 through 0{(3^^). This should be compared with the 
corresponding figure: 1.1 • 10^, in the LW computation. Since these figures by no means 
represent our computational limits, a further extension of our calculations is feasible. We are 
confident that our results are correct also because, in each space dimension d = 1, 2, 3, by 
a single procedure, we produce numbers in agreement with all expansion coefficients already 
available for A^ = 0,1,2,3 and oo (spherical model) Our codes were run on an IBM 



Rise 6000/530 power station with 32 Mbyte memory capacity and 1.5 Gbyte of disk storage. 
Typical cpu times were extremely modest and the RAM is far from being saturated. For 
reasons of space neither a detailed discussion of the main steps of this computation can fit 
here, nor we can display the extensive formulas giving the closed form structure of the HT 
series coefficients as functions of N. Therefore, as an example of our results, we shall only 
report here the HT series in the = 3 case (classical Heisenberg model) on the sc and 
the bcc lattices respectively, where we have contributed from five to eight new coefficients 
beyond those given in Refs. [p|,|9|,[T3| 



The susceptibility HT series for the sc lattice is: 

, 10^2 244^3 230^4 37612 . 864788^6 19773464^7 

y*=(/?) = 1 + 2/5 H 3^ H 3^ H 3^ H 3^ H 3^ H 3^+ 

^ ^' ' ^ 3 ^ 45 ^ 27 ^ 2835 ' 42525 ^ 637875 ^ 

89686514 ^„ 25478812^. 140348301868 477383158731608 
1913625 360855 1326142125 3016973334375 

426768736125964 12 28560817226680664 ,3 775988604270248 ^4 
1810184000625 ^ ^ 81458280028125 ^ ^ 1491909890625 ^ ^ 

16004552656617124832 ,5 354950851980427607594 
20771861407171875 ^ 311577921107578125 ^ ^ 



1464128352813955096312676 17 1068764655864454858376417828 
870237133653465703125 ^ ^ 430767381158465523046875 ^ 

259814093690188797550933157144 
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71076617891146811302734375 
The second correlation moment HT series for the sc lattice is: 
scfo^ a 02 964^3 2192^4 57116^5 8340368^6 33324872^7 1263947744^8 

73478278372 13325285538064^1. 3434294378983784^,. 51819882101501984^12 
63149625 6630710625 1005657778125 9050920003125 

773005999283909656 ^3 19031243835736702816 ^4 225650609227937809568 ^5 
81458280028125 ^ 1221874200421875 ^ 8902226317359375 ^ 
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267912260258927725784384 le 171544906778131647970688684 
6543136343259140625 ^ ^ 2610711400960397109375 ^ ^ 

45158335170568649028207863344 18 19874973349328684680550746792 ^9 
430767381158465523046875 ^ ^ 119456500657389598828125 ^ 

The susceptibility HT series for the bcc lattice is: 

^ 1 + 8^ + ^ 1^^3 + 12904 119600 2784992 632918848 
X ^p; ^ 3^ 9 ^ 135 ^ 405 1701 ^ 18225 1913625 ^ 

4075984504^0 287925718448 64384719769312 186782368415874752^,, 
5740875 ^ 189448875 ^ 19892131875 ^ 27152760009375 ^ 

1186773786369487616 12 7528780320376815776 ^3 732954612970918048 ^4 
81458280028125 ^ ^ 244374840084375 ^ ^ 11278838773125 ^ ^ 

127972570589148818590048 15 807291210775528531339816 
934733763322734375 ^ ^ 2804201289968203125 ^ ^ 

4736622265468109492081181616 17 1639367056527449858924222363488 ig 
7832134202881191328125 ^ ^ 1292302143475396569140625 ^ ^ 

56693738330512585673 156861 101868^ 
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213229853673440433908203125 

The second correlation moment HT series for the bcc lattice is: 

8^ 128^2 784^3 67072^4 4081648^,5 167636864 . 911026304^,7 

n'^J'^ff) = -3-\ 3^ H 3^ H 3^ H B-" H 3'' H 3^+ 

P2 v/^; 3/^ g -^5 g5Q5 127575 ^ 273375 ^ 

16849951744^8 4153759481008^0 1065794492624896 1166772237247486528^1, 

3 H H /3 H /? + 

1913625 ^ 189448875 ^ 19892131875 ^ 9050920003125 ^ 

8314233519972990976 12 175799675893471696544 13 409170117445661176448 14 
27152760009375 ^ ^ 244374840084375 ^ ^ 244374840084375 ^ ^ 

3613059270200364483884384 15 173915360520409186670373376 ig 
934733763322734375 ^ ^ 19629409029777421875 ^ ^ 

3161 1058478034436738314658288 17 84383259864054068055963337861 12 ig 
1566426840576238265625 ^ ^ 184614591925056652734375 ^ ^ 
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21963940464232232523449671606261888 19 
213229853673440433908203125 ^ 

Let us now comment on our updated estimates for the critical temperatures and the 
critical exponents 7 and v in the N = 2,3,4... cases where our new series are significantly 
longer than those previously available. The main difficulty of the analysis here comes from 
the expected singular corrections (conffuent singularities) to the leading power law be- 
havior of thermodynamical quantities. For example, the susceptibility should be described, 
in the vicinity of its critical point f3c, by 

X{N; f3) ^ C{N){(3, - (1 + a^{N){f3, - + ... + a'{N){f3, - f3) + ...) 

with the universal (for each N) exponent 0{N) ^ 0.5 for small N |jl^ and ^(A^) = 
1 + 0{1/N) for large N [^]. The standard ratio and Fade approximant (FA) methods are 
insufficient to cope with this very unstable double exponential fitting problem. Therefore 
we have to resort also to the inhomogeneous differential approximants (DA) method |2^, a 



generalization of the FA method better suited to represent functions behaving like 0i (x) (x — 
Xo)~'^ + 4'2{x) near a singular point xq, where 0i(x) is a regular function of x and 02(x) may 
contain a (conffuent) singularity of strength smaller than 7. We have essentially followed 
the protocol of series analysis by DA's suggested in Ref. |^ which is unbiased for conffuent 
singularities. We have computed Pc and 7 from the susceptibility series and have used this 
estimate of (3c to bias the computation of u from the series for the square of the (second 
moment) correlation length The results are reported in Table I along with the previous 
estimates by other methods. In the sc lattice case our exponent estimates are consistent 



with the RG e— expansion results ||T4|,|15[, but they are slightly larger (by ^ 1%) than the 
(7— expansion results. In the bcc lattice case the estimates are perfectly compatible with the 
most recent seventh order [|I^ or sixth order [|T^ (7— expansion results. This is analogous to 



what is observed in the most accurate unbiased analyses of the = 1 case and suggests 
that the series for lattices with lower coordination number have a slower convergence and 
also that unbiased DA's might be unable to account completely for the conffuent singularities. 
For A^ > 3 no elaborate estimates of the exponents by the e— expansion method are available 
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and only very recently an extensive computation by the (sixth order) expansion method 
has been published [jl^. Unfortunately, no estimates of error for the exponents are given in 
Ref. |]T^ , but we can safely assume uncertainties of the order of 0.5% for moderate values 
of and possibly smaller for > 8. 

Analysing our sc series by the simplest biased PA method designed to account explic- 



itly for the confluent singularities, does not significantly alter our DA estimates. Therefore, 
in order to assess with a higher level of precision the influence of these confluent singularities 
and completely reconcile the series results with those from the RG, further work is required 
including the computation of even longer series and, as indicated by the experience with 
the A^ = 1 study of suitable continuous families of models for each universality class 



rT| , |28[| . The uncertainties we have quoted for our exponent estimates, generously allowing 
for the scatter of the results in the DA analysis, leave small differences between our central 
values for the sc lattice, and those from the fixed dimension RG. This suggests that the still 
insufficient length and/or the still incomplete account of the confluent singularities add to 
some of our estimates a systematic uncertainty twice as large as we have indicated. Even 
under this conservative proviso, we have significantly improved the precision of the values of 
the critical parameters from HT series and have not pointed out any serious inconsistency 
with the estimates either from RG or from stochastic simulations. 
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TABLES 





TABLE I. A summary of the estimates of critical parameters 


for various 


values of N 


N 


Method and Ref. 


13c 


7 


V 


2 


Exper. 1^ 






0.6705(6) 




HTE sc 


0.45420(6) 


1.328(6) 


0.679(3) 




HTE bcc 


0.320434(8) 


1.323(2) 


0.674(2) 




HTE fee 


0.2075(1) 


1.323(15) 


0.670(7) 




J 1 

R.G. g-exp. |l6| 




1.318(2) 


0.6715(15) 




R.G. e exp. 10 




1.315(7) 


0.671(5) 




MonteCarlo se Q 


0.45420(2) 


1.308(16) 


0.662(7) 




MonteCarlo sc 


0.4542(1) 


1.316(5) 


0.670(7) 


3 


HTE se 


0.69302(7) 


1.403(6) 


0.715(3) 




HTE bee 


0.48681(2) 


1.396(3) 


0.711(2) 




HTE fee 


0.3149(6) 


1.40(3) 


0.72(1) 




J 1 

R.G. g-exp. |l6| 




1.3926(20) 


0.7096(15) 




R.G. e exp. |||| 




1.39(1) 


0.710(7) 




MonteCarlo se ^ 


0.693035(37) 


1.3896(70) 


0.7036(23) 




1 — 1 

MonteCarlo bee |32| 


0.486798(12) 


1.385(10) 


0.7059(37) 


4 


HTE se 


0.93582(8) 


1.471(6) 


0.749(4) 




HTE bee 


0.65526(3) 


1.458(3) 


0.742(2) 




R.G. g-exp. ^ 




1.45(3) 


0.74(1) 




R.G. g-exp. 




1.449 


0.738 




1 — 1 

MonteCarlo se [g3| 


0.9360(1) 


1.477(18) 


0.7479(90) 


6 


HTE se 


1.42838(9) 


1.577(6) 


0.801(4) 




xllll; Dec 


u.yyouo^^D ) 


1.564(3) 


0.795(2) 




R.G. g-exp. 




1.556 


0.790 


8 


HTE se 


1.9262(3) 


1.656(6) 


0.840(4) 




HTE bee 


1.33976(8) 


1.641(3) 


0.832(2) 




R.G. g-exp. 




1.637 


0.830 
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HTE se 


2.4267(3) 


1.712(6) 


0.867(4) 




HTE bee 


1.6850(2) 


1.696(3) 


0.859(2) 




R.G. g-exp. 




1.697 


0.859 
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